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I. INTRODUCTION 

Recent years have seen a rapid increase in the interest 
in finite quantum spin systems, also known as molecular 
magnets or nanomagnets Q&HSISIiSllSIIil- 
Molecular magnets typically consist of clusters of inter- 
acting spins that are magnetically isolated from the other 
clusters in the molecular solid by nonmagnetic ligands. 
Formally, molecular magnets are materials in which the 
ground state has nonzero total spin. Here we general- 
ize this definition to include all systems of largely iso- 
lated clusters of interacting quantum spins. These ma- 
terials are interesting both as simple model systems for 
the study of quantum magnetism and because they have 
possible applications as nanoscale computer memory el- 
ements [jj U- Many realizations of finite spin clusters 
with various ionic spins, ground state spins and geome- 
tries have been reported in the literature; some recent 
examples with S=l/2 ions are given in Table |I] 

Theoretical results for the properties of finite S=l/2 
quantum spin systems have appeared in several re- 
cent references, primarily in the context of experimen- 
tal studies of specific materials. Dimcr results arc 
reported in several studies of the S=l/2 spin dimer 
VO(HPO4)-0.5H2O; see for example Johnson et al. [lTI| . 
Tennant et al. [l3 and Koo et al. 0] . Theoretical prop- 
erties of S=l/2 spin trimers have similarly been given in 
studies of candidate trimer materials; see for example 
Refs.0HH[H[Ll|. 

Rather few general theoretical results have been re- 
ported for S=l/2 spin tetramers, since the results are 
more complicated and there are many more indepen- 
dent geometries and sets of superexchanges. Specific 
cases of tetramers are considered by Procissi et al. 
[T^| (S=l/2 square tetramer), Gros et al. [2(| and 
Jensen et al. [2lJ (an unsymmetric S=l/2 tetrahe- 



dral model of Cu 2 Te 2 5 (Bri_ a; Cl a; )), Kortz et al. 
(unsymmetric tetramer model of KyNafCu/tKafHaO^fq- 
AsWg033)2]-5.5H 2 0), and Ciftja 10] (symmetric trimer 
with apical spin). More general reviews of quantum 
spin systems have been published by Kahn [22| (ther- 
modynamics) and Whangbo et al. [23 (local origins of 
magnetism, thermodynamics properties, and materials). 
Studies of the dynamics of Heisenberg spin clusters using 
a quasiclassical formalism have been reported in a series 
of papers by Ameduri, Efremov and Klemm 0, US Ell- 
Waldmann [23 has carried out calcuations of the inelas- 
tic neutron structure factor for cyclic Heisenberg spin 
clusters which are quite similar to the results presented 
here. 

This increased level of interest in molecular magnets 
motivates more detailed theoretical investigation of the 
properties of finite quantum spin systems. For simple 
theoretical models such as the Heisenberg model, clus- 
ters that consist of only a few interacting magnetic ions 
can be treated analytically, and closed-form results can 
be obtained for many physical observables. One espe- 
cially interesting quantity is the inelastic neutron scatter- 
ing structure factor, which is required for the interpreta- 
tion of inelastic neutron scattering experiments. Inelastic 
neutron scattering is very well suited to the investigation 
of magnetic interactions at interatomic scales, since the 
measured structure factor is sensitive to the local geom- 
etry and interactions of the magnetic ions. As this work 
is intended in part to facilitate future neutron scattering 
studies, the evaluation of this structure factor is one of 
our principal concerns. 

In this paper we specialize to magnets that are clusters 
of S=l/2 ions with isotropic Heisenberg interactions, and 
give analytic results for the properties of dimer, trimer 
and tetramer clusters with various geometries. After the 
introduction, in Sec. II we define the Heisenberg model 



2 



T Jr\ 13 ±j J_j 1 . OUHlfc; cA.dlllJJlt'b 


f^T email ^% — 1 /V riiiavi+nm cmn cvrctoTnc 
UL Ollld'll O — J- / £t LJ LlctllL LIII1 OjJlll OVfaLCIllB. 








ojjiii oysLem 


Prnnnrl ^ 
VjrlUUIlO DLdLc Otot 


IV CIS. 


VO(HPO 4 )-0.5H 2 O 


dimer 







CU3(02Cl6H23)6 ■ 1.2CeHl2 


symmetric trimer 


1/2 


[16, 17] 


Nag [Cu 3 Na 3 (H 2 0) 9 (a-AsWg0 33 )a] -26H 2 


symmetric trimer 


1/2 


m 


[Cu 3 (cpse)3(H 2 0) 3 ] ■ 8.5H 2 


symmetric trimer 


1/2 


[31] 


(CN 3 H6)4Na 2 [H4V608(P04)4((OCH 2 ) 3 CCH 2 OH) 2 ]-14H 2 isosceles trimer 


1/2 


14] 


Na 6 [H 4 V 6 8 (P04)4((OCH 2 ) 3 CCH 2 OH) 2 ]-18H 2 


general trimer 


1/2 




K 6 [Vi5As 6 042(H 2 0)]-8H 2 


symmetric trimer + capping hexamers 


1/2 


[32, 33, 34, 


NaCuAs0 4 


linear tetramer 





[28, 29J~ 


(NHEt 3 )[Vi 2 As 8 4 o(H 2 0)]-H 2 0. 


rectangular tetramer + capping tetramers 





[37] 


K 7 Na[Cu 4 K 2 (H 2 0) 6 (a-AsWg0 33 ) 2 ]-5.5H 2 


distorted tetramer 


1 


w 



and the observables we evaluate in this work. These in- 
clude the standard thermodynamic quantities for mag- 
netic materials (partition function, specific heat and mag- 
netic susceptibility) , as well as the inelastic neutron scat- 
tering structure factors. In Section III we evaluate these 
quantities for specific spin clusters, which are the spin 
dimer, symmetric, isosceles and general spin trimers, and 
three cases of spin tetramers (tetrahedral, rectangular, 
and alternating linear). We also tabulate all the energy 
eigenvalues and eigenvectors for each spin system. Sec- 
tion III ends with an application, which is a numerical 
study of powder inelastic neutron scattering amplitudes 
in NaCuAs04 ,28]; the results appear to support the 
identification of this material with the alternating linear 
tetramer model and agree well with the data of Nagler 
et al. 29] . Finally, Section IV discusses several mate- 
rials which may be candidates for future experimental 
studies, as well as some interesting directions for future 
theoretical research. 

For reference purposes our principal results for the spin 
systems considered here are given in a series of tables 
at the end of the paper. These results are the spec- 
trum of energy eigenvalues and eigenvectors (Table HH. 
allowed inelastic neutron scattering transitions between 
these states fTable lTTlj) . and the specific heats f Tabic ITV|) 
and susceptibilities (Table [Vj. 

Most previous theoretical studies of molecular magnets 
in the literature have specialized to individual materials 
and their associated model Hamiltonians. Our results are 
intended to be sufficiently general so that they should be 
useful for the interpretation of data on many candidate 
molecular magnets. 



nets, is defined by the Hamiltonian 



n = 



(i) 



where the superexchange constants {Jy } are positive for 
antiferromagnetic interactions and negative for ferromag- 
netic ones, and Si is the quantum spin operator for a 
spin- 1/2 ion at site i. 

Since this is a rotationally invariant Hamiltonian in 
spin space, the total spin Stot is a good quantum num- 
ber. For the specific cases of dimer, trimer and tetramer 
clusters of S=l/2 ions that we consider here, the energy 
eigenstates have the total spin decompositions given be- 
low. 



1/2® 1/2 = 18 0, (2) 



1/2 ® 1/2® 1/2 = 3/2 ffil/2 2 , (3) 



1/2® 1/2® 1/2® 1/2 = 2© 1 3 ©0 2 . (4) 

Each Stot multiplct contains 2Stot + 1 magnetic states, 
which are degenerate given an isotropic magnetic Hamil- 
tonian such as the Heisenberg form of Eq.JTJ. 



B. Expressions for Observables 



II. THE MODEL AND OBSERVABLES 

A. The Heisenberg Magnet 

The nearest- neighbor Heisenberg magnet, which we 
shall assume as our standard model for molecular mag- 



The energy eigenstates and eigenvalues may be found 
by diagonalizing the magnetic Hamiltonian on a conve- 
nient basis. (In practice we will employ the usual set of 
z-polarized magnetic basis states.) Several physically in- 
teresting quantities may be computed directly from the 
energy eigenvalues; in this work these are the partition 
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function, specific heat and magnetic susceptibility, which 
are given by 



Z = ^e-^=^(2S tot + l)e-^, (5) 

i=l Ei 



C = k B (3 



2 d 2 HZ) 
d[3 2 ' 



(6) 



and 



N 



2^. e -/3E, 



J(.9Mb) 2 §E( 2S ^ 



l)(S tot + l)S tot e-' 3E ' . (7) 



In these central formulas the sum i = 1 . . . JV is over all 
N independent energy eigenstates (including magnetic 
substates), the sum J2e ^ s over energy levels only, M z — 
mgfiB where m — S% ot /H is the integral or half-integral 
magnetic quantum number, and g is the electron g-factor. 

In addition to these bulk quantities, we also give results 
for inelastic neutron scattering intensities. In "spin-only" 
magnetic neutron scattering at zero temperature, the dif- 
ferential cross section for the inelastic scattering of an in- 
cident neutron from a magnetic system in an initial state 
l^i), with momentum transfer hq and energy transfer hw, 
is proportional to the neutron scattering structure factor 
tensor 



/CO 7 i 
?L ^ e *(^-^)+-*(* i |St(x j ,t)S a (x i ,0)|* i ) . 

(8) 

The site sums in Eq.lJSJ run over all magnetic ions in 
one unit cell, and a, b are the spatial indices of the spin 
operators. 

For transitions between discrete energy levels, the time 
integral gives a trivial delta function <5(Ef — E; — hui) in the 
energy transfer, so it is useful to specialize to an "exclu- 
sive structure factor" for the excitation of states within 
a specific magnetic multiplet (generically ^/(A/))) from 
the given initial state 

Sj£\?) = E PW^/CV)) (*f(*f)\Va\*i) , (9) 



where the vector V a (q) is a sum of spin operators over 
all magnetic ions in a unit cell, 



(10) 



This exclusive structure factor is related to the exclu- 
sive differential inelastic neutron scattering cross section 
by 

Ml = ( 7 r ) 2 |(<5 Qb - ms { £>{<tm?)\ 2 (ii) 

where 7 w —1.913 is the neutron gyromagnetic ratio, 
r = ha/m e c is the classical electron radius, k and k' are 
the magnitudes of the initial and final neutron wavevec- 
tors, and F(q) is the ionic form factor. (This relation is 
abstracted from Eq.(7.61) of Ref . [3fT| . specialized to an 
exclusive process.) 

For a rotationally invariant magnetic interaction and 
an St t = initial state (as is often encountered in 
T=0 inelastic scattering from an antiferromagnet), only 

Stot = 1 final states are excited, and S b ^\q) cx 5 a b. 
In this case we may define a scalar neutron scattering 
structure factor S(q) by 



S ( b f J ) (q) = S ab S(q) 



(12) 



The result for Sf£ (q) is more complicated for neutron 
scattering from a magnetic (Stot > 0) initial state. If we 
assume an isotropic magnetic Hamiltonian and a spheri- 
cal basis for the spin operators S a , the tensor Sf£ % \q) is 
diagonal but is not cx S a b] it instead has entries that are 
proportional to a universal function of q times a product 
of Clebsch-Gordon coefficients, since 

(*f(\ f )\ V a |^(A0) = (Sf\f\la Sj Aj) V^(q) , (13) 

where V^ 1 ' (q) is the reduced matrix element for the tran- 
sition — > |\&/). Here we simplify the presentation by 
quoting the unpolarized result (Sj£ ((f)), obtained by 
summing over final and averaging over initial polariza- 
tions. This unpolarized (Sj£ z \q)) is cx S a b, so it suffices 
to give the function S(q); 

(Sl f : ) (q))=5 ab S(q) = 



2Si 



1 ^ 



(*i(A,)|V 6 t |\£7(A / )) (* / (A / )|K|*i(Ai)) . 

If desired, the general results for polarized scattering can 
be recovered by reintroducing the appropriate Clebsch- 
Gordon coefficients of Eq.JTSJl in Eq.@. 

The results given above apply to neutron scattering 
from single crystals. To interpret neutron experiments 
on powder samples, we require an orientation average of 
the unpolarized single-crystal neutron scattering struc- 
ture factor. We define this powder average by 



S(q) 



dQ r ~, 



4tt 



1 S(q, 



(15) 
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III. RESULTS FOR SPECIFIC CASES 
A. Spin Dimer 
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FIG. 1: The geometry and energy levels of a Heisenberg spin 
dimer. 



The "minimal" spin cluster model is the S = 1/2 spin 
dimer (Fig^J, which consists of a single pair of S = 1/2 
spins interacting through the Heisenberg Hamiltonian, 



Ji = J Si • S2 . 



(16) 



Since this is an isotropic magnetic Hamiltonian, the total 
spin is a good quantum number, and from the Clebsch- 
Gordon series 1/2 £g) 1/2 = 1 © we expect the spectrum 
to consist of an Stot = 1 triplet and an Stot = singlet. 
In a i-diagonal basis 



in: 
in: 
in: 
in: 



(17) 



the Hamiltonian matrix is 
1/4 



n 



-1/4 1/2 
1/2 -1/4 



1/4 



(18) 



Diagonalizing this Hamiltonian matrix gives the energy 
eigenvalues and eigenvectors, 



Ei = 
E = 



(19) 
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0)) = 






-1)) = 


III) 



ITT)) 



l*o) 



V2 



(ITI)-IIT)) 



(20) 



(21) 



The specific heat and magnetic susceptibility for the 
dimer are especially simple, since there is only a single 
excited level. The results are (in a dimensionless form) 



c_ 

k B 



3(/3J) 2 



-pi 



(l + 3e-^Y 



and 



X 



Cwb) 2 /J 



2/3J 



(l + 3e-^ J ) 



(22) 



(23) 



(24) 



These results are summarized in Tables IIVI and as 
are the corresponding results we find for the other spin 
systems considered in this paper. Plots of the dimension- 
less specific heat and susceptibility of the spin dimer are 
shown in Figs l2l3l 




FIG. 2: The magnetic contribution to the specific heat of a 
spin dimer, Eq. l23H (dimensionless units). 

One may confirm that this specific heat formula gives 
the correct entropy for a dimer of S=l/2 ions, 



S= /°°C7^ = fc B 21n(2) . 
Jo P 



(25) 
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FIG. 3: The magnetic susceptibility of a spin dimer, Ea. (|24ll 
(dimensionless units). 



FIG. 4: A fit of the dimer susceptibility formula of Eq. l24H to 
the measured susceptibility of VO(HPO4)-0.5H2O. A defect 
term was also included. 



The corresponding result for a general spin system is 

S = k B HAf/Afo) , (26) 

where N is the dimensionality of the full Hilbert space 
and No is the degeneracy of the ground state; for the 
S=l/2 dimer, N = 2 2 and N = 1. 

As an example of the application of the dimer suscep- 
tibility of Eq. 1)24|) (known as the Bleaney-Bowers formula 
|38|). in FigQ]we show a fit to the susceptibility of the 
spin dimer VO(HPO 4 )0.5H 2 O HJ. (The molar suscep- 
tibility shown is related to the single dimer susceptibility 
of Ea.lJ23|l by Xmoiar = ■ X-) The parameters of the 

fit are g = 2.05 and J = 7.76 meV (consistent with the 
results of inelastic neutron scattering [l2j ) . A 1/T defect 
contribution was also included in the fit. 

Finally we evaluate the inelastic neutron scattering 
intensities, which are given by the structure factors of 
Eas. l|14ll5|l . (A complete set of inelastic neutron scat- 
tering transitions for all the spin systems we consider in 
this work is given in Table IIHI typically we will only 
evaluate the structure factors for the ground state of the 
antifcrromagnetic system.) We evaluate Ea. (|14|l for the 
dimer using the energy eigenvectors |\I/i(m)) and \^o) °f 
Eqs.HOEIl)- This gives 

S(q) = ±(l-cos(q-a)) (27) 

where a = x\ — x^ — x\i is a spatial vector that coincides 
with the dimer. Evidently there should be no excitation 
of the dimer spin-triplet state when the neutron momen- 
tum transfer q is perpendicular to the dimer axis a. 

In scattering from powder samples one measures the 
powder average S(q) of the structure factor, defined by 
Ea. p5|l . For the dimer this is 

S(q) = J ^ ~ (1 - cos(g • a)) = \ (l - j (qa)) (28) 



1 -0 I — i — i — i — i — ■ — i — ■ — i — ■ — i — ■ — i — ■ — i — ■ — i — ■ — r 
0.8 - 




2 4 6 8 10 12 14 16 18 20 



qa 

FIG. 5: The powder average, unpolarized neutron structure 
factor S(q) for a spin dimer with pointlike magnetic ions. 

where jo (x) — sin(x) /i is a spherical Bessel function. 
This result is shown in Fig0for pointlike magnetic ions 
(F(q) — 1). The location of the first maximum, at 
q fa 4.493 a -1 , provides a convenient estimate of the sep- 
aration between the interacting ions in the dimer. Of 
course in real materials the incorporation of ionic form 
factors will reduce the location of this maximum. 

Experimental studies of real magnetic materials typi- 
cally proceed by establishing the approximate magnetic 
parameters of a model Hamiltonian through a fit to 
the susceptibility. Given a model Hamiltonian, one can 
predict the inelastic neutron scattering structure factor, 
which is then compared to experiment. (Ideally this is 
done on single crystal samples, but frequently only pow- 
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der samples are available.) Unlike the bulk susceptibility, 
the inelastic neutron scattering structure factor allows a 
sensitive and microscopic test of the assumed magnetic 
Hamiltonian, since it is determined by the relative posi- 
tions of the interacting magnetic ions. The spin-dimcr 
material VO(DPC>4)-0.5D2O provides a recent illustra- 
tion of the use of inelastic neutron scattering in identi- 
fying magnetic interaction pathways; the susceptibility 
data of Johnson et al. 1 1 1 1 | was well known to give an 
excellent fit to the dimer formula Ea. (|24ll . however the 
separation of the interacting V-V pair inferred from in- 
elastic neutron scattering data ^3 using Eq. I|28[) showed 
that the interacting V-V pair had been misidentified in 
the literature. 



B. Trimers 

We will consider the most general case of a spin trimer 
with Hciscnberg magnetic interactions. It is useful to 
present the results as special cases with decreasing sym- 
metry, since the formulas arc simpler in the more sym- 
metric cases, and examples of both symmetric and isosce- 
les trimers are known in the literature. 



1. Symmetric Trimer 

The completely symmetric, equilateral trimer has 
equal magnetic couplings and bond lengths between all 
three pairs of spins. The Hamiltonian for this model is 



T~L — J I Si • S2 + S2 • S3 + S3 • Si 



(29) 



Since this Hamiltonian is invariant under any permuta- 
tion of the three spin labels, it has a discrete S3 symmetry 
in addition to the magnetic rotational symmetry. In an 
S z -diagonal basis 

TTT 
Hi 
TIT 
ITT 
Til 
III 
III 



(30) 



the Hamiltonian matrix is 



3/4 



-1/4 1/2 
1/2 -1/4 



1/2 
1/2 



1/2 1/2 -1/4 



-1/4 
1/2 
1/2 



1/2 
4/4 
1/2 



1/2 
1/2 
-1/4 



3/4 



(31) 

This matrix is block diagonal within subspaces of defi- 
nite S 2 tot, as expected for a rotationally invariant Hamil- 
tonian. The energy levels of the symmetric trimer are 
shown in FigEJi. For the J > (antiferromagnetic) case 
the ground state is a quadruplet (the two Stot = 1/2 mul- 
tiplets are degenerate), and there is an energy gap of |J 
to the S( t = 3/2 excited state. Representative symmet- 
ric trimer energy eigenstates (those with maximum S^ tot) 
are given in Table ITT1 Since the two Stot = 1/2 levels are 
degenerate, there is no unique ground state for this sys- 
tem; we use the Jacobi \p) = \{12)a) and |A) = |(12)s) 
three-body basis states of definite (12)-exchange symme- 
try as our two independent basis vectors. 

We may determine the specific heat and magnetic sus- 
ceptibility of the symmetric trimer from these energy lev- 
els, using Eas. l|6l7fl . The results are 



C_ 

k B 



W\2 



and 



X _ 1 (l + 5e-|/») 
(g» B ) 2 /J 4 P (l + e-^ J ) 



(32) 



(33) 



It is notable that the integral of this specific heat gives 
an entropy of 



S = /°°C'^ = A :i3 ln(2) 
Jo P 



(34) 



which is only half as large as the entropy of the dimer, 
despite the larger trimer Hilbert space, Af = 2 3 = 8. The 
lower entropy is due to the fourfold degenerate ground 
state of this highly frustrated system; 

S = k B HAf/Afo) = k B ln(2 3 /4) = k B ln(2) . (35) 

The susceptibility of the symmetric trimer, Eg 1331 
agrees with Eq.2 of Veit et al. 01 (after specializing 
to a single g-factor and a change of variables). This re- 
sult is shown in Fig[7| Note that x(T) diverges as we 
approach zero temperature, since the ground state has 
nonzero spin. This divergence is present independent of 
whether the intrinsic spin-spin coupling J is antiferro- 
magnetic (as we normally assume) or ferromagnetic, since 
both cases have ground states of nonzero spin. A more 
detailed comparison of the susceptibility suffices to dis- 
tinguish these; see the inset of Fig[7| which shows %T ver- 
sus T for both cases. At high temperatures the spin-spin 
coupling J is unimportant, and both results approach the 
same Curie's law limit. 

Next we consider the neutron scattering structure fac- 
tors for the symmetric trimer. Since this system has two 
degenerate Stot = 1/2 ground states and a single St D t = 
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FIG. 7: The magnetic susceptibility of a symmetric trimer. 
The inset shows yT versus T for ferromagnetic (dashed) and 
antiferromagnetic (solid) couplings. 



3/2 excitation, there are two distinct exclusive inelastic 
neutron structure factors but only a single transition en- 
ergy, Ei — E = | J. We have chosen \p) and |A) basis 
states for our orthogonal S tot = 1/2 eigenstates, and will 
give neutron structure factors for each of these. The same 
structure factors follow for the isosceles trimer, although 
in that case the two Stot = 1/2 states are nondegenerate. 

The structure factors for excitation of the Stot = 3/2 
| a) level (using Ea. lfHjl ) are given by 



1 



- (l + -cos(g-»i 2 )- 



• cos(g-xi 3 ) 



iq-x 23 )) • (37) 



These results may be understood in terms of the dif- 
ferent natures of the \p) and |A) initial states. In the 
\p) ground state the (12)-dimer is in a pure S(i2) = 
state, which must be excited to S(i 2 ) = 1 to couple to 
the |er) excited state. The \p) — > \<r) excitation problem 
is thus identical to the dimer problem, to within an over- 
all constant. It follows that S^ p ^ a '(q) is proportional to 
the dimer structure factor of Ea. H27|) . In contrast, in the 
|A) initial state the (12)-dimer is pure S(i 2 ) = 1 and the 
(23)- and (31)-dimers have amplitudes to be in both spin 
and 1, so there are contributions to S^^^ due to the 
excitation of each of the three dimer subsystems. 

As ( S , (^ cr ) and 5( A_+cr ' differ considerably for moderate 
qa it will certainly be possible to distinguish between \p) 
and | A) states from their single crystal structure factors. 
The powder averages however are identical, and cannot 
be distinguished experimentally; 



1 



S^°\q) = S^°\q) = U\-j {qa)) 



(38) 



S^(q) = -{l-cos(q-x 12 )) , 



(36) 



These powder structure factors are identical because 
of the identical dimer lengths, n 2 = r 2 3 = ?*3i = a, so 
the powder average of each cosine in Eas. (|36l37|l gives the 
same ja{qa) Bessel function. The dependence (1— jo{qa)) 
follows from the requirement that 5(0) = 0. 

As we shall discuss in the next section, an isosceles 
trimer would be a more favorable system for the identi- 
fication of \p) and | A) initial states in inelastic neutron 
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scattering; these levels are nondegenerate in the isosceles 
system, and the \p) — ► \a) and |A) — * \a) powder average 
structure factors are no longer equal, due to the different 
leg lengths. 



2. Isosceles Trimer 

The isosceles spin trimer, Fig. [SJd, has two equal mag- 
netic interactions and bond lengths. The Hamiltonian is 
given by 



TL = J( Si • S2 + a(S2 • S3 + S3 • Si ; 



(39) 



To find the energy eigenvalues of this Hamiltonian it 
suffices to consider the S z tot = +1/2 sector, since the 
Stot = 1/2 and 3/2 multiplets both have S ztot = +1/2 
members. The remaining symmetry of this problem sug- 
gests that we use the three {|S tot , +1/2)} energy eigen- 
states of the symmetric trimer as our basis, 

H3/2, +1/2)) = yfl (| m ) + 1 TIT ) + 1 ITT )) 
|A(l/2,+l/2)) = yj\ CI TXT > + IITT) - 2|TU» 
| P (i/2,+i/2)) = ^T (ITIT)-IITT)) 

The Hamiltonian is necessarily diagonal on this basis, 
since these three basis states have different values of the 
conserved quantities Stot and (12)-exchange symmetry. 
The result is 



n 



1 + 2a 



1 - 4a 



(41) 



The two S tot = \ levels are split as a result of the re- 
duced symmetry of the isosceles trimer; the full S3 sym- 
metry of the symmetric trimer has been reduced to S2 
((12)-exchange symmetry), and since S2 is Abelian no 
degeneracies follow from this symmetry. 

The specific heat and susceptibility of the isosceles 
trimer, which follow from the energy levels of Ea. l41|) 
and the formulas Eas. l|t)l7|) . are given in Tables ITVl and 
IVl The susceptibility agrees with the earlier result of 
Veit et al. [4(|. Note that one recovers the symmetric 
trimer result in the limit a = 1. 

We have confirmed by numerical integration of the 
rather complicated isosceles trimer specific heat formula 
given in Tabic ITVI that the entropy of the isosceles trimer 
satisfies 



S 
k B 



d£ 



21n(2) 
ln(2) 



a ^ 1 
a = 1 



(42) 



as expected from Eq. Ij26(l for an eight-dimensional Hilbcrt 
space which has a fourfold-degenerate ground state for 
a = 1, and a twofold-degenerate ground state otherwise. 



Although the magnetic contribution to the specific 
heat is usually masked by much larger phonon contribu- 
tions, we note in passing that one may separate magnetic 
contributions experimentally by subtracting the specific 
heats in zero and nonzero magnetic fields. This approach 
was used recently by Luban et al. 0] to study an S=l/2 
V 4+ vanadium trimer, (CNaHe^^pHUVeCVPO^ 
((OCH 2 ) 3 CCH 2 OH) 2 ]-14H 2 (their material 1), which 
appears to be an accurate realization of the isosceles 
Heisenberg trimer. 

As we found for the symmetric trimer, the susceptibil- 
ity of the isosceles trimer also diverges as T approaches 
zero, since the system has a magnetized ground state. 
The rate of divergence with T can again be used to distin- 
guish between ferromagnetic and antifcrromagnetic cou- 
plings (which have St„t = 3/2 and Stot = 1/2 ground 
states respectively), as shown in Fig|7|for the symmetric 
trimer. This behavior is evident in the susceptibility of 
material 1 of Luban et al. [l4|; in Fig. 3 of this reference 
one can see that %T for this material clearly follows the 
lower trimer curve, confirming that it is accurately de- 
scribed by the antiferromagnetic isosceles trimer model 
(with an Sf t = 1/2 ground state). 

There are three inelastic transitions excited by neutron 
scattering from an isosceles spin trimer, \p) — > |er), |A) — > 
\a) and \p) — ► |A). The first two were considered in the 
discussion of the symmetric trimer, and the results for 
the isosceles trimer are identical (except that the AE 
values of the transitions differ). The \p) — * |A) transition 
was not considered previously because these states are 
degenerate in the symmetric trimer. The result we find 
for the structure factor of this transition is 

S^(q) = ±(l-co S (q-x 12 )) . (43) 

This has the same form as the dimer and \p) — > |er) struc- 
ture factors because it also involves the excitation of the 
S(i2) = (12)-dimer to an S(i2) = 1 state. It can evi- 
dently be distinguished from the \p) — > \a) transition by 
the overall intensity, but not by the functional depen- 
dence on q. 

To illustrate these single crystal structure factors, in 
Fig|5]we show the two ground state structure factors for 
the previousl y c ited V6 isosceles trimer (material 1 of 
Luban et al. [14()- The parameters are a = 3.22 A and 
b = 3.364 A. We show the predictions of Eas. l|37l43|l for 
in-plane scattering, with momentum transfer q = n/a. 
Since this material has two strong bonds and a weak 
dimer (a ~ 9) 0], M should be the ground state, and 
the I A) — > \p) and |A) — > \cr) transitions shown in the 
figure should both be observable (These are expected at 
4.2 meV and 7.0 meV respectively, given the parameters 
of Luban et al.) The very different angular distributions 
predicted for the scattered neutrons show that it should 
be straightforward to distinguish between these transi- 
tions in an inelastic neutron scattering experiment, given 
a single crystal of this or a similar trimer material. 

The powder average eliminates much of the difference 
between these neutron scattering transitions, although it 
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- — ■ C 




FIG. 8: The unpolarized structure factors {S(q}} (propor- 
tional to the angular scattering intensities) predicted for in- 
elastic neutron scattering from the [A) ground state of a single 
crystal of an isosceles trimer material (see text). 



still should be possible to distinguish them experimen- 
tally. On carrying out the powder average we find 



|(i -jo(qa)) 

|(l + |io(go)-|jo(«6)) (44) 

i(l-joM) • 



In the symmetric limit b/a = 1 these transitions are 
proportional to the same function, 1 — jo(qa); at best it 
may be possible to distinguish the \p) — > |A) transition 
from the others through their relative intensities. How- 
ever for significantly different leg lengths the |A) — > \a) 
powder average structure factor may differ enough from 
the 1 — jo(qa) of the \p) — > \a) and \p) — » |A) transitions 
to distinguish them. As an example, in FigOwe show the 
powder structure factors of Eq. 144JI for the three transi- 
tions, for an elongated triangle with b/a = 2. (These 
results are independent of the magnetic coupling ratio 
a.) As there is considerable variation in form and mag- 
nitude between these powder structure factors, it should 
be possible to distinguish them experimentally in similar 
isosceles trimer materials. If more than one transition is 
clearly observed, it may also be useful to compare struc- 
ture factor ratios, to eliminate the effect of ionic form 
factors. 



3. General Trimer 

The general trimer of FiglU has three different mag- 
netic couplings and ion pair separations, and is described 
by the Hamiltonian 



H = J( Si • S 2 + «23 §2 • S3 + "31 S3 • Si 



(45) 
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FIG. 9: The powder average inelastic neutron structure factor 
S(q) for the three allowed transitions of an isosceles trimer, 
with b/a = 2. 



This system is also discussed by Qiu et al. JT3 in the 
context of La4Cu3MoOi2, which they model as a two 
dimensional coupled array of S=l/2 trimers. 

We may again determine all the trimer energy eigen- 
values by specializing to the S ztot = +1/2 sector and 
using the symmetric trimer basis of Eq. I|4(J|) , which gives 
the Hamiltonian matrix 



H = - J 

4 



1 



1 — 2a s V3ctd 
V3ad —3 



(46) 



where a s — a^i + a.23 and ad = CX31 — «23- The \a) basis 
states again must be energy eigenstates, since they are 
the only Stot = 3/2 states in the Hilbert space. They 
have energies of 



Es = 



1 



(47) 



The Stot = 1/2 basis states \p) and |A) mix in this prob- 
lem, since the general trimer Hamiltonian with «23 7^ c*3i 
breaks (12)-exchange symmetry. The resulting energies 
are 



Ei 



(l + a 8 ±v/(2-a a ) 2 +3^) 



(48) 



The specific heat and susceptibility of the general 
trimer follow from the energy levels of Eas. l|47l48|l and 
the formulas Eas. H6l7|) . The resulting expressions are 
given m Tables Hvl and M One may confirm recovery 
of the isosceles and symmetric trimer results as special 
cases of these results. We have also confirmed by numer- 
ical integration that the rather lengthly general trimer 
specific heat formula given in Table IIVI leads to an en- 
tropy of S = fee 2 ln(2), provided that at least one of the 
parameters «23 and 031 differs from unity. 
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The neutron scattering structure factors for the general 
trimer involve coherent superpositions of the previously 
derived \p) and |A) excitation functions, since the energy 
eigenstates are superpositions of these basis states. The 
Stot = 1/2 energy eigenstates of Eg . 1)46(1 are explicitly 



•sin(0)|A) + cos(0)|p) 



and 



\* h2 } = +co S (d)\\)+sm(9)\p) 
where the mixing angle satisfies 
tan(6>) = f 

i + yiT^2 



(49) 



(50) 



(51) 



with x — V3ad/(2 — a s ). The Stot = 3/2 energy eigen- 
state is, as for all the trimers we have considered, 



1*3 



(52) 



The structure factor for the transition from the Stot = 
1/2 state |*i x ) to the Stot = 3/2 state is given by 



^2 ( C31 



C23 — 2 C12) 



C23) (53) 



where dj = cos(<f • scy), C 2 — cos(20) and S 2 = sin(20). 

The structure factor for the second transition, * 1 9 — > 

2 ' z 

*3 , follows from Eg. (|53l> on changing the overall signs of 
the C 2 and S 2 terms. The third transition, between the 
two S tot = 1/2 states, has the structure factor 



S 



'an 



1 



1 



(l--(Cl2 + C 2 3 + C 3 l)) 



+ -C4 (C31 +C23 



2 C12) H — j=. S4 (C31 — C23) 
V3 



(54) 



where the new quantities are C4 = cos(4#) and £4 = 
sin(40). One may confirm that the previously de- 
rived symmetric and isosceles trimer structure factors of 
Eas. l|36l37l43|) follow from these general trimer results in 
the limit 6 -> Q. 

The powder averages of these general trimer unpolar- 
ized structure factors may also be evaluated; the result 

for the transition *i , — * *3 is 

2 2 

^i^!),. iA (1 + 2C 2 ) . , > 



(1 - C 2 + v 7 ^) . , , (1 - Ca - v 7 ^) . , . 
5 Jo (qr 2 3) 5 Jo (qr 3 i ) 



The powder average results for the two remaining tran- 
sitions can be obtained from Eq. H55f) by simple substitu- 



te* 1 



tions. To obtain S 5 simply change the overall 

signs of C-y and So in Ea. H55(l . and to obtain S , 
divide Eq. 1)55)1 by a factor of two and replace C 2 and 1S2 
by C4 and S4 respectively. 

These results will be useful for the interpretation of 
neutron scattering data on real materials. One ex- 
ample of a candidate general trimer is the V§ mate- 
rial 2 of Luban et al. [lj, Na 6 [H 4 V 6 8 (P04)4((OCH2)3 
CCH20H)2]T8H20. This compound has three distinct 
V-V separations between the S=l/2 V 4+ ions within each 
vanadium trimer, 3.212 A, 3.252 A and 3.322 A. 



C. Tetramers 

We will consider three S = 1/2 tetramer spin clusters of 
decreasing symmetry, the regular tetrahedron, the rect- 
angular tetramer, and the linear (dimer-pair) tetramer. 
Our definitions for the magnetic couplings and geometry 
of these systems are shown in FigllOl As with the dimer 
and trimer systems we will give results for the partition 
function, specific heat, magnetic susceptibility and neu- 
tron inelastic scattering structure factors, the latter for 
both single crystal and powder average cases. 



1. Tetrahedron 

This system has four S = 1/2 ions at the vertices of 
a regular tetrahedron, with Heisenberg interactions of 
strength J between each pair of ions (see FigllOfa). The 
Hamiltonian of this system is given by 



Ti = J S ' Si ■ Sj 



(56) 



i<3 



The invariance of this Hamiltonian under permutation 
of any site labels implies an S4 symmetry, in addition 
to the spin rotation symmetry SU(2). Since the group 
S4 is non-Abelian and has irreducible representations of 
dimensionality 1, 2 and 3, we anticipate that one may 
find twofold and threefold degeneracies in the spectrum 
of tetrahedron energy eigenstates. We will see that this 
is indeed the case. 

As with the dimer and symmetric trimer we may de- 
termine the energy eigenvalues of this system by simply 
squaring the total spin operator Stot = Yh=i which 
gives for this case 



E s t . 



IK 



Stot (St. 



(55) 
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FIG. 10: Energy level diagrams for (a) the tetrahedron, (b) the rectangular tetramer, and (c) the linear (dimer-pair) tetramer. 



The Clebsch- Gordon series of Eq.Q implies that these 
Stot = 1 and Stot = energy levels are respectively three- 
fold and twofold degenerate. 

Given these energy levels, the specific heat and suscep- 
tibility of the tetrahedron may then be determined using 
Eas. l|t)l7|) . with the results 



C ? ,,(1 + 50 

k B 2 KH ' f-< . 9„ 



-2/3J 



- 10e- 3 ^ J ) 



(l + |e-/3.T + | e -3/3J) 2 



and 



(3Ms)7J 



3/3Je 



-03 _ 



(1 



-2/3J 



1 



-/3J i 5 p -3f33) 



(58) 



(59) 



These quantities are shown in Fig^Jand Figll2lrespec- 
tively. The specific heat of the tetrahedron gives an en- 
tropy of S = ks 31n(2), as expected for a 16-dimensional 
Hilbert space and a doubly-degenerate ground state. 

Note that the susceptibility is rather similar to that of 
the spin dimer, since the tetrahedron also has an Stot — 
ground state and a gap of J to the magnetized Stot = 1 
excited states. (The fact that the ground state is twofold 
degenerate does not affect this result, since both are 
Stot = states and neither makes a contribution to the 
susceptibility.) 

Determination of the energy eigenvectors requires di- 
agonalization of the Hamiltonian on a specific basis. Op- 
erating on our | (12) (34)) dimer basis of Eas. l|A.3IA.5|) 
with the tetrahedron Hamiltonian, Ea. (|56|l . we find that 
the Hamiltonian matrix is already fully diagonal; each of 
these basis states is an energy eigenvector of the tetrahe- 
dron Hamiltonian. 

In our discussion of neutron scattering structure fac- 
tors of the tetrahedron and the other spin tetramers con- 
sidered in this paper, we will specialize to Stot = initial 




k B T/J 

FIG. 11: Magnetic contribution to the specific heat of a reg- 
ular tetrahedron. 



states. Structure factors for Stot > initial states, which 
are of interest for systems with magnetized ground states 
and at finite temperatures, and can be derived using sim- 
ilar methods. 

The tetrahedron has two degenerate ground states, 
which we take to be l^o.i) = \pp) and 1^0,2) = |cc)o- 
The three degenerate Stot = 1 excited states, which 
can be reached from the S tot = levels using inelas- 
tic neutron scattering, are taken to be = \(pa)s}, 
|*i,2) = \(pct)a) and |*i, 3 ) = \(aa) s )i- The choice of 
this specific set of initial and final states is rather arbi- 
trary; in a real material we would expect a spontaneous 
distortion of the lattice, which would select nearly degen- 
erate energy eigenstates that need not be these specific 
basis states. However these will suffice to illustrate the 
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FIG. 12: Susceptibility of a regular tetrahedron. 



neutron scattering structure factors expected for nearly 
tetrahedral systems. 

The single crystal structure factors for all of these tran- 
sitions may be read directly from Eas. ijA.lllA.18fl . For 



example, the transition |\&o,i) 
matrix element of Ea. (|A.ll|) ; 
definition in Ea. (|14f> . we find 

5 (*o,i-*i,i)(^) = !/ 2 



— > l^i.i) is specified by the 
using the structure factor 



- (C12 — C13 + C14 + C23 — C24 + C34) / 4 



(60) 



where as before dj = cos(q ■ 2%). This characteris- 
tic angular distribution and its five partner distributions 
could be used in an inelastic neutron scattering exper- 
iment from a single crystal sample to characterize the 
spin states of the individual Stot = and Stot = 1 lev- 
els. (Note however that one specific transition, |\Po.i) — ¥ 
I '5 1.3), has a zero matrix element.) 

The powder average structure factors for a tetrahe- 
dron are much less characteristic. Since there is only a 
single ion pair separation, each cosine in the single crys- 
tal structure factors such as Ea. (|6U[) powder-averages to 
the same factor of jo(qa). This gives a powder struc- 
ture factor that is proportional to 1 — joiqa) for each 
transition, just as we found for the dimer and symmet- 
ric tetramer; only the overall coefficients distinguish the 
different transitions. These results are 



gOo, i-»*i,i 


(<?) = 
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(1 


-joiqa)) 


C(*0,l->*1,2 


(<?) = 
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(1 


-joiqa)) 


g^O, 1-»*1,3 


(<?) = 









g(*0, 2-**l,l 
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(1 


-jo{qa)) 


5(*o,a-»*i,2 


(<?) = 
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(1 


-joiqa)) 


g(*0, 2-**l, 3 


iq) = 


I 


(1 


-joiqa)) 



(61) 



A generalization of the tetrahedron problem in which 
the Hamiltonian has couplings of strength aJ between 
ions in different dimers may also be of interest. This 
generalized Hamiltonian is 

H{a) = jf(Si-S 2 + S 3 -S 4 ) 



(62) 



a ( Si • S3 + Si • S4 + S2 • S3 + S2 • S4 



The dimer-pair basis of Eas. (|A.3IA.5|) is also diagonal 
under this Hamiltonian, with the eigenvalues given be- 
low. (We have added an Stot subscript to all these state 
vectors for clarity.) 



(63) 
(64) 
(65) 
(66) 
(67) 
(68) 



- -ha J \aa) 2 
1 

a 

2 

1 



aah 



H(a)\a<r) S{ot =2 
W(a)k<7)i 
Hia)\(pa) s )i 
Hia)\ipa) A ) 1 

Hia)\aa) = - 2a) \aa) 
H(a)\pp) 



-5 \(P a )s)i 

~\ \iP0)A)l 



\PP)( 



Since the energy eigenvectors of this generalized prob- 
lem are exactly the basis states we used for the tetra- 
hedron, the neutron scattering structure factors for the 
Stot = to Stot = 1 transitions are unchanged. In 
this system however all these levels are nondegenerate, so 
unlike the pure tetrahedron problem one encounters no 
structure factor ambiguities due to an arbitrary choice 
between degenerate basis states. 



2. Rectangular Tetramer 

The rectangular tetramer, shown in FigllUb. has (12) 
and (34) dimers of interaction strength J coupled by in- 
teractions of strength aJ between ion pairs (13) and (24). 
The Hamiltonian is 



H = J( Si • S 2 + S 3 • S 4 + a(Si • S 3 + S 2 • S 4 ) 



(69) 



This Hamiltonian is already diagonal on the Stot = 
1,2 dimer-pair basis states of Eas. (|A.3IA.4|) ; the energy 
eigenvalues are 



E'2 

Ei, 3 

El,2 

Ei 1 



(1 
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(1 
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(1 
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J . 



(70) 
(71) 
(72) 
(73) 
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The 2x2 Hamiltonian matrix in the Stot 
spanned by Ea. HA.5fl is 

-3 



1 - 2a 
— Voa 



subspace 



(74) 



which has the eigenvalues 



E, 



o,{>} 



1 1 / 

- — -a =f \/l — a- 



(75) 



The specific heat and susceptibility of the rectangu- 
lar tetramer may be evaluated using these energy levels 
and the general formulas of Eas. (|6l7|) . The susceptibil- 
ity is given in Table [V] Although the specific heat is 
straightforward to evaluate, the resulting expression is 
too lengthly to tabulate here. 

The neutron scattering structure factors of the rect- 
angular tetramer are especially interesting because the 
ground state is a linear combination of the dimer-pair ba- 
sis states; this mixing leads to coupling-dependent struc- 
ture factors. The S tot = ground state of the rectan- 
gular tetramer is a linear superposition of unexcited and 
doubly-excited dimer pairs, 



|*o,i> = -sin(0 o )|<7o-)o +cos(9 )\pp) 



(76) 



where the mixing angle 6$ between these basis states sat- 
isfies 



tan(#o) 



V3a/2 



1 - a/2 + Vl - a + 



(77) 



The matrix elements of the neutron scattering tran- 
sition operator V a of Eg. (|1U|) between the ground state 
|\&o. i) and the three Stot = 1 final states |^i,i...3) may 
then be determined using the results of Eas. l|A.lllA.18|) . 
Specializing to S z tot — +1 final states for illustration, 
these matrix elements are 

l,l(*l,„|V+|*0,l> = 



(Co+Sq/V3) (f 

n = 3 , 



/2 + /S-/4), 

■% (A + h-h- fi), 



(78) 



where Co and So are cos(#o) & n d sin(#o) respectively. 

On converting these matrix elements to structure fac- 
tors using Eas. 1)12114(1 we find different functional forms 
for each transition for both the single crystal and powder 
average results, as a result of the different weight factors 
and the three distinct ion pair separations. The powder 
average, unpolarized structure factors are 

C*0,l->*1, 



'(?) 



(Cq+Sq/%/3)" 



(1 



3o{qa)- j (qb) +jo(qc)), n = 1 
jo (qa) + jo (qb) - j (qc)) , n = 2 
jo(qa)~ jo (qb) - jo {qc)) , n = 3 

(79) 
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FIG. 13: Powder average, unpolarized structure factors S(q) 
for the excitation of the three S tot = 1 excited states of the 
rectangular spin tetramer from the Stot = ground state 
l^o.i)- The interdimer coupling strength is a = 0.3, and the 
side length ratio is b/a — 1.5. The structure factor for the 
transition to the third state is scaled up by a factor of 10 for 
visibility. 
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FIG. 14: Variation of the S(q) structure factors for excitation 
of the lowest Stot = 1 state of the rectangular spin tetramer 
with side length ratio b/a. This illustrates the use of powder 
structure factors in establishing the internal geometry of spin 
clusters. The dimer coupling is a — 0.3, however this only 
affects the overall normalization. 



where c = v 'a 2 + b 2 . In Fig El we show these struc- 
ture factors for a case with moderate interdimer coupling 
(q = 0.3) for a rectangle of side ratio b/a — 1.5. The 
transition to the highest S tot — 1 excited state ^1,3) is 
much weaker than the other two, and so is multiplied by 
a factor of 10 in the figure for visibility. 
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Note that the excitation of the highest Stot — 1 state 
1*1,3), which is a doubly-excited dimer (\aa)), is only 
possible because the ground state has an O(a) excited 
component, in addition to the dominant "bare" \pp) ba- 
sis state. The weakness of the 1*1,3) signal is because the 
structure factor is proportional to the nonleading ground 
state amplitude squared, so that it is 0(a 2 ). The obser- 
vation of similar "non-valence state" transitions which 
are forbidden at 0(a°) should allow direct experimental 
tests of the "interaction" terms in quantum spin Hamil- 
tonians such as this one. The strength alone is a sen- 
sitive measure of a, and the spatial modulation of S(q) 
is clearly different from the 0(a°) transitions that dom- 
inate the structure factors of the lower-lying Stot = 1 
states. 

The detailed q-dependence of the structure factors can 
be used as a "fingerprint" to test whether a given struc- 
ture is indeed the magnetically active system. As an 
example, in Fig^Jwe show that the detailed form of the 
structure factor to the first Stot = 1 state shows signif- 
icant variation with the ratio b/a. (Values of b/a = 1, 
2 and 4 are shown.) This type of dependence could be 
used to establish the geometry of a magnetic subsystem, 
or to check powder neutron scattering results against the 
geometry of a proposed spin system. 



Linear Tetramer 



with eigenvalues 



E 



Mil 



-I -a ± -y/l + a 2 ) J 



(84) 



and 



E 



o,C 2 } 



\ + \a±Jl-\a + \ a > J (85) 



respectively. 

The susceptibility of the linear tetramer, which follows 
from these energy levels and Eq.{7|), is given in Table [V] 
As with the rectangular tetramer, the expression we find 
for the specific heat is too lengthly to tabulate here. 

The neutron scattering structure factors from the lin- 
ear tetramer S to t = ground state 1*0,1) to the three 
Stot = 1 excited states |*i,i...3) may be calculated us- 
ing the same techniques we applied to the rectangular 
tetramer. The results are somewhat more complicated, 
since two of the S tot = 1 linear tetramer states are ro- 
tated between the dimer-pair basis states, in addition to 
the ground state basis rotation we found for the rectan- 
gular tetramer. The energy eigenvectors in these sectors 
are the superpositions of dimer-pair basis states given in 
Table ITT1 with mixing angles that satisfy 



The linear tetramer consists of two dimers with inter- 
nal magnetic couplings of strength J, with a single inter- 
dimer coupling of strength aJ between the two adjacent 
end spins (see Fig llUfl . The term "linear" refers only to 
the pattern of magnetic couplings; the actual spatial ge- 
ometry of our linear tetramer is not assumed to be a 
straight line. The "Clemson tetramer" NaCuAsOi [2^] 
is a recent example of a possible "linear tetramer" that 
does not have a true collinear dimer geometry. 

The linear tetramer Hamiltonian matrix is also rela- 
tively simple in the dimer pair basis of Eas. (|A.3IA.5|l . 
The single S tot = 2 basis state \(aa)}2 is diagonal, with 
energy 



E? 



1 1 



a) J 



(80) 



,2 4 

The Stot = 1 basis state \{p<j)a) is also diagonal, with 
energy 



Ei, 2 = -It-tq|.i 



(81) 



The remaining Stot = 1 and Stot = two-dimensional 
Hamiltonian matrices are 



Hs tot =i - - J 



1-ia 



and 



H Stot =o - 2 J 



#a -3 



(82) 



(83) 



tan(0i) = 



l + Vl + a 2 



tan( 



V^a/4 



1 - a/4 + y/l -a/2 + a 2 /4 



(86) 
(87) 



This more complicated basis mixing pattern introduces 
a new feature, which is that the functional forms of the 
structure factors for the two mixed Stot = 1 states, |*i,i) 
and 1*1.3), depend on the dimer coupling a. (In the rect- 
angular tetramer system discussed previously we found 
that changing the dimer coupling a only changed the 
overall normalization of the structure factors, not their 
detailed q dependence.) 

We will give explicit results for the first transition, 
|*o,i) - * l*i,i), and then simply quote the results for 
the two remaining final states. The matrix element of 
the neutron scattering transition operator V a is 



11 (*l,l|V+ |*0,l) = C U (ft - h) + C23 (/2 - / 3 ) 

where 

SqSi 



C14 
C23 



±- 



2V2 



V6 



(89) 



and Co,i and 5o,i are respectively cos and sin of the linear 
tetramer basis state mixing angles, which were defined in 
Eqs. jH3S3 . 
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FIG. 15: Powder average, unpolarized structure factors S(q) 
for the excitation of the three Stot = 1 states of the lin- 
ear tetramer, including the ionic form factor, with magnetic 
coupling ratio a — 0.4 and ion positions taken from the 
NaCuAsC>4 structure |2S|. The structure factor for the tran- 
sition to the third state is scaled up by a factor of 10 for 
visibility. 



FIG. 16: Contours of equal intensity of the powder average, 
unpolarized structure factors of the linear tetramer, with pa- 
rameters appropriate to NaCuAs04, as in Fig ll5l The con- 
tours are normalized to the peak of the first transition, and 
intervals of 0.2 in intensity are shown for the two lower states. 
The much weaker third transition is displayed with intervals 
of 0.01. 



The resulting unpolarized powder structure factor for 
this transition, using Eas. (|14ll5fl . is 



'(^ui 1 _ io(gi"u)) + c 2 2 3 (i - io(gr 2 3)) 



+2ci 4 c 23 (jo(<?T 2 ) - kiOTxs)) 



(90) 



The transition from the ground state to the second 
linear tetramer Stot — 1 excited state is given by a ma- 
trix element we encountered previously in the rectangular 
tetramer problem, except for a change in spatial geom- 
etry. The result for the unpolarized powder structure 
factor with completely general ion positions is 



(Co + S /Vsf 



+jo(<7?"23) - jo(^24) - jo(qru)) 



(91) 



The structure factor for the third S tot = 1 state ^1,3 
may be found from the structure factor of Ea. H90|) 
with the simple substitutions (G\ — > Si), (Si — > —C\). 

We will illustrate the predicted structure factors for 
the linear tetramer assuming parameters appropriate for 



the candidate material NaCuAs04 28, 29] . The ion pair 
separations are r%2 — = 3.641 A, ri3 = r 2 4 = 3.863 A 
ri4 = 6.814 A and r 2 3 = 3.151 A. The copper posi- 
tions are consistent with planarity, although our results 
do not require this assumption. There are indications 
of three Stot = 1 levels in this material from a recent 
inelastic neutron scattering experiment, with energies 
of approximately 9, 11 and 18 meV 29]; in the linear 
tetramer model this suggests parameters of J « 10 meV 
and a ~ 0.4, which we will assume here. (The structure 
factors only depend on a.) We also incorporated a sim- 
ple Cu 2+ ionic form factor, F(q) = 1/(1 + q 2 /^) 3 with 
go = 8.0 A -1 , which agrees with the online ILL Cu 2+ 
form factor to < 0.5 % over the range of q shown. 
Our results are shown in Fig^J characteristic features 
include the displaced relative maxima of the intensities 
of the two lower states, and the much weaker transition 
to the highest Stot = 1 state. The results for the two 
lower states are rather insensitive to a. The overall scale 
of the 1^1,3) structure factor however is quite sensitive 
to a, and scales approximately as a 2 . A measurement of 
the relative strength of these transitions would provide a 
useful determination if a, which could be compared with 
the value extracted from the energy levels. (In principle 
the susceptibility could also be used to determine a, but 
we have found that it has a rather weak a dependence in 
this system.) 

In Fig^lwe show these results in a contour plot, ap- 
proximately as would be observed in a neutron scatter- 
ing experiment. (Our results should be multiplied by 
the energy-dependent factor fe'/fc of Eq.JTTJ for a direct 
comparison with experiment.) To generate this plot we 
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have convolved a Gaussian energy resolution function, 
exp(-(E - Ei) 2 /2a%) with a E = 0.5 meV, with the 
structure factors to the three Stot — 1 states. The in- 
tensities are shown relative to the maximum excitation 
intensity of the transition to the lowest Stot = 1 state, 
l^o, l) - * l^i, i) ■ Note the characteristic strong peak in 
intensity of the second transition, I^Vi) — ¥ 1^1.2), near 
1.1 A -1 . Comparison of these general features with the 
data of Nagler et al. suggests that the linear tetramer 
model does indeed give a realistic description of neutron 
scattering from NaCuAsO-4. 



IV. FUTURE APPLICATIONS 

In the previous section we presented results for 
bulk thermodynamic and magnetic properties of vari- 
ous dimer, trimer and tetramer molecular magnets with 
S = 1/2 ions. We also derived the inelastic neutron scat- 
tering structure factors for these systems; inelastic neu- 
tron scattering is very useful as a local probe of magnetic 
interactions at the atomic scale. These results clearly 
have many possible applications to real materials. In 
this section we discuss some examples of materials that 
are thought to be realizations of S = 1/2 dimer, trimer 
and tetramer molecular magnets, and describe how our 
results could be useful in future experimental investiga- 
tions. We also discuss possible extensions of this work, 
which will be useful in interpreting experimental data on 
these and related magnetic materials under more general 
circumstances. 

The S = 1/2 spin dimer is the simplest of all spin clus- 
ters. It provides a textbook case for studies of finite spin 
systems more generally, since physical observables for the 
dimer can often be derived as closed form analytic expres- 
sions. This is the case for the specific heat, susceptibility 
and neutron scattering structure factors presented here. 
Vanadyl hydrogen phosphate, VO(HPO4)-0.5H 2 O, is a 
well known example of an S = 1/2 spin dimer material 
[TH IT2I H^ l , and some of the results tabulated here have 
already been used in interpreting data on this material. 
In particular, the magnetic susceptibility was originally 
used to determine the exchange constant J, and inelas- 
tic neutron scattering was used to test the simple dimer 
model and establish which pair of V 4+ ions forms the 
dimer |T^ . This experiment was a dramatic success for 
inelastic neutron scattering, as the previously assumed 
V-V dimer pair was shown to have been misidentified. 

Many examples of S=l/2 ion trimers have been re- 
ported in the literature. These systems are interesting in 
that the ground state is (ideally) degenerate, and must 
exhibit fcrromagnetism. (Stot > for any energy eigen- 
state of an isotropic magnetic Hamiltonian with half- 
integer ion spins and an odd number of ions.) Heisen- 
berg trimers with antiferromagnet pair interactions are 
also of interest because they are the simplest isotropic 
spin systems which experience frustration. One example 
of an S=l/2 trimer is Cu 3 (0 2 Ci 6 H 23 ) 6 -1.2C 6 H 12 OGll, 



which has an equilateral Cu 2+ triangle with a Cu 2+ 
- Cu 2+ separation of 3.131 A. Recent Electron Para- 
magnetic Resonance (EPR) measurements show that 
the ground state of this material consists of a twofold- 
degenerate Stot = 1/2 level; this is in accord with ex- 
pectations for a general isotropic trimer antiferromagnet 
with S = 1/2 ions, but not with the perfect equilateral 
(symmetric) case, in which the ground state is a quar- 
tet of two degenerate Stot = 1/2 levels. The gap to the 
Stot = 3/2 excited level is estimated from both EPR and 
susceptibility data to be 28 meV 0, . There are in- 
dications from the EPR studies that this fourfold degen- 
eracy has been lifted by additional, nonisotropic interac- 
tions 01 ■ Investigation of the level structure and struc- 
ture factors in this apparently symmetric trimer material 
would be a very interesting exercise for a future inelastic 
neutron scattering experiment, especially if large single 
crystals are available. 

Another recent example of an S=l/2 trimer is the 
"Na 9 -2" material of Kortz et al. [3, Na9[Cu 3 Na 3 (H 2 0)9 
(a-AsWg0 33 ) 2 ]-26H 2 0, which contains an equilateral 
Cug + trimer with a susceptibility consistent with equal 
Heisenberg interactions of J ~ 0.35 meV. Inelastic neu- 
tron scattering from a powder sample of this material 
should show the Stot = 3/2 excited level, with a struc- 
ture factor proportional to (1 — jo(qa)). With a single 
crystal sample it might be possible to separate the tran- 
sitions from the two (nearly?) degenerate Stot = 1/2 
levels to the S to t = 3/2 excited level. 

The two "V 6 " materials of Luban et al. 14], (CN 3 H 6 ) 4 
Na 2 [H 4 V 6 8 (P0 4 )4 ((OCH 2 ) 3 CCH 2 OH) 2 ] -14H 2 and 
Na 6 [H 4 V 6 8 (P04)4 ((OCH 2 ) 3 CCH 2 OH) 2 ]-18H 2 0, con- 
tain pairs of (presumably weakly coupled) V 3 spin 
trimers that are respectively isosceles and general tri- 
angular systems. The isosceles Vg material was used as 
an example of single-crystal inelastic neutron scattering 
structure factors in this paper. 

The "Vi 5 " material K 6 [Vi 5 As 6 42 (H 2 0)]-8H 2 is an 
example of an S=l/2 trimer in a more complicated mag- 
netic geometry. This material has a frustrated V 3 trian- 
gle sandwiched between two nonplanar antifcrromagnetic 
V6 hexagons. The low-temperature magnetic properties 
are dominated by the V 3 triangle; other magnetic in- 
teractions become important at elevated temperatures 
IH S S |M E3. In addition to distinguishing 
between direct vanadium-vanadium and superexchange 
pathways involving the upper and lower hexagons, neu- 
tron scattering was used to probe the magnetic structure, 
finding two nearly degenerate S{ t = 1/2 ground states 
(with 0.035 meV splitting) and an S to t = 3/2 excited 
state [35J . The availability of large single crystals of this 
material suggests that it might be an interesting candi- 
date for inelastic neutron scattering studies. 

Clearly, the analytical expressions we have presented 
here for spin-trimer thermodynamic properties and in- 
elastic neutron scattering amplitudes have wide potential 
application, and should be useful in particular for inter- 
preting the results of future inelastic neutron scattering 
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experiments on spin-trimer molecular magnets. 

Examples of tetramer systems with S = 1/2 ions in- 
clude sodium copper arsenate, NaCuAsC>4 H^,|2i|, which 
we used as an illustration of the evaluation of inelas- 
tic neutron scattering structure factors in the previous 
section. The neutron scattering data of Nagler et al. 
|29j | supports a model of this material as an open-chain 
tetramer, with antiferromagnctic Hcisenberg bonds of al- 
ternating strength. Transitions from the Stot — ground 
state to all three Stot = 1 triplet excited states have been 
observed on a powder sample j2^| , and the relative inten- 
sities and the g-dependence of the powder average struc- 
ture factors appear to be approximately consistent with 
our predictions for the open-chain model. As we have 
given detailed analytic predictions for the neutron struc- 
ture factor for these transitions, a comparison with data 
from a high-statistics experiment on a larger powder sam- 
ple should be straightforward. 

Additional examples of S = 1/2 tetramers are found in 
vanadium materials containing the [Vi 2 Ass04o(H 2 0)] 4 ~ 
cluster anion. Basler et al. |37j ] have recently reported 
studies of the magnetic properties of three such materials, 
Na 4 [Vi 2 As 8 4 o(H 2 0)]-23H 2 0, Na 4 [Vi 2 As 8 04o(D 2 0)] 
•16.5D 2 and (NHEt 3 )[Vi 2 As 8 4 o(H 2 0)]-H 2 0. These 
materials have three stacked V4 tetramers, but are 
mixed- valent (Vg + V^ 5 ). The middle tetramer domi- 
nates the magnetic properties. This tetramer is anti- 
ferromagnetic and close to square, with exchange 
constants of ~ 1.5 meV (inferred from energy levels 
established by EPR and inelastic neutron scattering). 
The Basler et al. study is a very nice illustration of the 
combined use of bulk magnetic properties and inelastic 
neutron scattering to characterize magnetic materials, 
as we advocate in this work. Additional studies of this 
already well characterized material might involve an in- 
elastic neutron scattering study of a single crystal, which 
could be used to test the detailed orientation dependence 
expected for the structure factor for each of the observed 
magnetic transitions to excited states, given their fitted 
magnetic Hamiltonian. Since this Hamiltonian includes 
anisotropies, a study using polarized neutrons could 
provide additional useful information. 

There are several interesting questions which were not 
considered in detail in this paper that would be appro- 
priate for future research on finite spin clusters. Con- 
sideration of higher ionic spin is one obvious general- 
ization of this work. Several examples of uncompen- 
sated molecular magnets (which have ground states with 
nonzero spin) may be found in relatively simple higher- 
spin materials. One example is the first cobalt molecu- 
lar magnet H3, Co 4 (NC5H4H 2 CO)4(CH30H) 4 Cl4. This 
material consists of four S = 3/2 Co 2+ ions and four 
ligand-related oxygen atoms situated on the corners of 
a cube, with a ferromagnetic Stot — 6 ground state. 
The magnetic exchange constants have been estimated 
from fits to magnetization curves, and are at the meV 
scale [43|. The magnetic Co 2+ ions in this material 
form a tetramer with important tetrahedral and dimer 



magnetic interactions |43j, which would be an impor- 
tant case for future neutron scattering studies, espe- 
cially with large single crystals. Another example of a 
molecular magnet with higher ion spin is the chromium 
magnet [Cr 4 S(0 2 CCH 3 ) 8 (H 2 0) 4 ](N03) 2 H 2 0, which has 
four ferromagnetically coupled S = 3/2 Cr 3+ ions ar- 
ranged in a nearly regular tetrahedron, with an Stot — 6 
ground state. Furukawa et al. JeJ determined the ex- 
change constant for this material from the susceptibility, 
and predict a gap to the first Stot = 5 excited state of 
ca. 15 meV. Observation of this excitation using inelas- 
tic neutron scattering should be a straightforward exer- 
cise, and the intensities should agree well with theoretical 
expectations of the Heisenberg model, since this model 
gives a reasonably good description of the bulk magnetic 
properties. 

Extension of this work to mixed-valent spin clus- 
ters would also be interesting, since many examples of 
these are known, including systems with magnetized 
ground states. One example is the "Mn4" material 
p^n 4 3 Cl 4 (0 2 CEt)3(py)3]2-2C 6 Hi4 studied by Hill et al. 
[44|. which has a mixed-valent spin tetramer consisting 
of a triangle of S = 2 Mn 3+ ions with an apical S = 3/2 
Mn 4+ , and a ground-state spin of Stot = 9/2. Although 
the interest in this material as a molecular magnet is 
largely due to the weak coupling between pairs of Mn 4 
clusters, the magnetic Hamiltonian within a single Mn 4 
cluster could be tested by powder inelastic neutron scat- 
tering. 

Another interesting theme for future studies is the ef- 
fect of finite temperatures on inelastic neutron scattering; 
although increasing temperatures are usually associated 
with weaker inelastic transitions, the finite Hilbert space 
of a spin cluster implies that magnetic transitions will 
weaken in a simple, known manner according to their 
Boltzman factors, and will approach finite limits at high 
temperatures (provided that the magnetic Hamiltonian 
remains valid). Since moderate temperatures (on the 
scale of the magnetic excitations) will significantly pop- 
ulate excited levels, it may also be possible to observe 
inelastic transitions from excited levels that are inacces- 
sible at low temperatures. 

Another important topic which we briefly alluded to 
in the text is the issue of magnetic interactions between 
spin clusters; these interactions will broaden the discrete 
levels assumed here into bands, which will be observable 
if the intercluster interactions are sufficiently large. 

Finally, the generalization of our results to non- 
Heisenberg interactions, and the determination of these 
interaction parameters through polarized inelastic neu- 
tron scattering experiments, would be an especially in- 
teresting and important extension of the work presented 
here. 
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SUMMARY AND CONCLUSIONS 



and 



In this paper we have evaluated several thermody- 
namic and neutron scattering observables that charac- 
terize the magnetic behavior of finite quantum spin sys- 
tems. After an introduction that gives results applicable 
to the general case, we specialized to clusters of S=l/2 
ions with a Heisenberg interaction between ion pairs. We 
considered dimer, trimer and tetramer systems with var- 
ious magnetic interaction strengths, and evaluated the 
magnetic specific heat, the susceptibility, and the inelas- 
tic neutron scattering structure factor for these systems. 
The structure factor was derived both for single crys- 
tals and for the powder average case. Our results for 
the neutron scattering structure factor show that accu- 
rate intensity measurements of inelastic neutron scatter- 
ing cross sections from a powder can be useful in estab- 
lishing the spatial geometry of an assumed set of inter- 
acting magnetic ions. The linear spin-tetramer candidate 
NaCuAsCU was considered as an example, and we found 
that the observed inelastic powder pattern for excitation 
of the two lowest Stot = 1 excited levels is indeed con- 
sistent with the predictions of the linear tetramer model. 
We also considered inelastic neutron scattering from sin- 
gle crystals, and found dramatic angular dependence that 
could be used in future experiments as sensitive tests of 
the assumed magnetic Hamiltonian. We concluded with 
a discussion of specific materials that might be studied 
using our results, and suggested future extensions of our 
work to more general systems. 
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These are combined as Clebsch- Gordon series to form 
tetramer basis states of definite total spin and symme- 
try, which are \(p<j) s ), S tot = 0,2; \(aa) A ), S tot = 1; 
\(p<?)s,a)i Stot = 1; \pp), S to t = 0. In the interest of clar- 
ity we will occasionally specify the total spin of one of 
these basis state with a subscript; thus \aa}o refers to 
the \(aa)s) state with S to t = 0. 

Using these states as basis vectors reduces the 16- 
dimensional full tetramer Hilbert space to 1-, 2- and 3- 
dimensional subspaces, which are spanned by the basis 
sets 



2) 



{ I Stot = 1) } = 



{ I Stot = 0) } 



\cra) 2 

\0*r)s) 
A(pvU) 

|cr<T> 
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(A.3) 



(A.4) 



(A.5) 



Thus symmetry arguments alone determine the eigenvec- 
tors for one level, and the eigenvectors for the remaining 
levels involve at most 2x2 and 3x3 diagonalizations. As 
we shall see, for the three tetramer models we consider 
here we actually encounter at most 2x2 diagonalization 
problems using this basis. 

These basis states are also convenient for determining 
neutron scattering structure factors, since they have rela- 
tively simple matrix elements of the spin transition oper- 
ator V a of Eq. I|1U|) . The complete set of matrix elements 
of V a (spherical components) between single (12)-dimer 
basis states, with fi = e ' , is 



APPENDIX: Tetramer basis states and matrix 
elements 

Since there is a natural separation of the rectangle 
and linear tetramer systems into dimer components, it is 
useful to introduce a | (12) (34)) dimer basis to represent 
tetramer energy eigenvectors. The dimer basis states are 



(p\V a \p) = 
(a(m) | V a \p) = &a, 



h-h 



(p | V a \cr(m)) = 



2 

fi-h 



(A.6) 
(A.7) 

(A.8) 



(a(m')\ V a \a(m)) = -5 m ,, m+a . (A.9) 



\P) 



1 

71 



(m>-in>) (a.i) 



These dimer results may be combined to give the com- 
plete set of matrix elements of V a between tetramer basis 
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states, which is all that we require to determine all neu- 
tron scattering structure factors for all the spin tetramer 
problems we consider. These tetramer matrix elements 
(with explicit Stot or Stot, Sztot subscripts on the states 
where required for clarity) are 



(pp\V a \pp)=Q 



l,m((P a )s I V a \pp) = 



S a ,: 



l,m((p<j)A | V a \pp) = 



Sa,i 



2V2 



h — h — h + h 



(aa I V a \pp) = 



2^2 



(A.10) 
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(A.12) 
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2 (aa I V a \aa) = 
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The remaining matrix elements between pairs of 
Stot = 1 states and between Stot = 1 and Stot = 2 states, 
which were not required in this paper, may be evaluated 
similarly. 



(pp\V a \aa}=0 
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TABLE II: Energy eigenvectors and eigenvalues. 



Spin System Eigenvector |^s to t) (Sztot = Stot) Energy 



Dimer 


l*i> = 


= k(+i)) = |TT) 


Ei = 








|*o) = 


= \p)= Tf(ITl) - lit)) 


Eo = 


3 T 




Symmetric Trimer 


1*3) 

2 


= K+3/2)) = ITTT) 


E 3 = 

2 


2 J 

4 J 






1*1 2! 

2 


= |A(+i/2)) = ^(|UT) + |ITT)-2|TU)) 


Ei = 

2 


4 J 






1*1 i) 


= |p(+i/2)> = ^(|UT>-UTT» 








Isosceles Trimer 


l*3> 

2 


= k(+3/2)> 


E 3 = 

2 








1*1 2 ! 


= |A(+l/2)> 


El 

2 ' Z 








1*1 l) 

2 ,1/ 


1 / 1 1 /o\\ 

= IM+ 1 / 2 )) 


El ! 


= ~4 J 




General Trimer 


1*3) 

2 


= k(+3/2)> 


E 3 = 

2 


= £(i + a.)J 








= +cos(0)|A(+l/2)> + sin(0)|p(+l/2)> 


Ei 

2 > J 


= i(+V(2-« S ) 2 + 3^ - 


l-a s )J 






= - sin(0)|A(+l/2)) + cos(0)|p(+l/2)> 




= i(-V(2-a s ) 2 + 3^ - 


l-a s )J 


Tetrahedron 


1*2) = 


= ka) a = ITTTT) 


E2 


2 J 






1*1,3) 


= k^)i 


1 = 


1 T 






1*1,2) 


= |(pff)s> 










1*1,1) 


= (p°")a) 










1*0,2) 


= ctct)() 


Eo = 








1*0,1) 


= M 








Rectangular tetramer 


1*2) 


= \cra) 2 


E2 


(i + ia)J 






1*1,3) 


= ko-)i 


f 


= (±-±a)J 






1*1,2) 


= l(^)s> 


Ei, 2 


= (-! + !«) J 






1*1,1) 


= I(^M 


Ei,i 


= (-±-i«)J 






1*0,2) 


= +cos(6» ) \aa)o + sin (# )pp) 


Eo,2 




J 




1*0,1) 


= — sin(0o)|crcr)o + cos(0o)\pp) 


Eo,i 


= (-Vl-a + a 2 -i-ia; 


J 
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TABLE II: Energy eigenvectors and eigenvalues (cont). 



Spin System 



Linear Tetramer 



Eigenvector |* Stot ) (Sztot = Stot) 



1*2) = |ctct) 2 

|*i,3) = +cos(0i)|cr<7)i +sm(e 1 )\(pa)s) 
|*1,2) = \( P *)a) 

|*i,i) = -sin(0i)|cr(7>i +cos(6»i) I (pa) s ) 
|*o,2) = +cos(6> )|(t<t)o + sin(0o)|pp) 
|*o,i) = -sin(0o)|<Tcr)o + cos(0 o )|pp) 



Energy 



E 2 = (I + ia)J 



Ei, 3 = 
Ei, 2 
f — 

Eo,2 = 

Eo 1 = 



-i + |a)J 



3«)J 



■iVT+?-ia)J 



+ 1 Q+ 1 Q 2_ 1 



7«) J 



- V /l -5«+3« 2 -|-I«)J 
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TABLE III: Inelastic Neutron Scattering Transitions" 



System 




Transition 


AE 


Dimer 


I. 


|*o) |*i) 


J 


Symmetric Trimer 


I. 


1*3,1/ 


_> lilr o \ 

— 1^3) 


2 J 




II. 




-|*J> 




Isosceles Trimer 


I. 


l* 5 ,l) 




(l + |a)J 




II. 




-» |* 3 ) 


faJ 




III. 


1*1,!) 




(l-a)J 


General Trimer 


I. 


I*i,l) 


-|*f) 


(I(l+a s ) + I 




II. 


l* 3 ,2> 


-|*|) 


(i(l+a s )-i 




III. 


1*1,!) 







Tetrahedron I. 


l*o,i) — 


|*i,i) 


J 


TT 
11. 


l*o,i> — > 


|*1,2) 




III. 


l*o,i) — > 


|*i,3) 




IV. 


1*0,2) — » 


|*l,l) 




V . 


1 iTf \ 

1*0,2/ — » 


|Wl,2; 




V 1. 


1*0,2/ — » 


1*1,3; 




Vll. 


1*1,1/ — 


1*2/ 


O T 

zj 


VTTT 
V 111. 


1*1,2/ — > 


1*2/ 




IX. 


*1 •?) — » 


* 2/ 




Rectangular Tetramer I. 


1*0.1) — 


l*l,l) 


/iJ 


II. 


|*0.l) - 


|*1,2) 


(/1 + «) J 


III. 


l*0,l) - 


1*1,3) 


(/1 + 1) J 


IV. 


l*l,l)- 


|*0,2) 


/iJ 


V. 


l*l,l)- 


1*1,2) 


a J 


VI. 


l*l,l)- 


*i,s) 


J 


VII. 


l*l,l)- 


1*2) 


(1 + a) J 


VIII. 


1*1,2) - 


1*0,2) 


Cfi - ") J 


IX. 


1*1,2) - 


*i,s) 


(l-a)J 


X. 


1*1,2) - 


1*2) 


J 


XI. 


|*0,2) -> 


1*1,3) 


(l-/i)J 


XII. 


1*1,3) - 


1*2) 


Q J 


Linear Tetramer I. 


l*0,l) - 


1*1,1) 


(/2-i(/3-l))J 


II. 


l*0,l) - 


1*1,2) 


(/a + ±a) J 

(/2 + i(/3 + l))J 
(/2 + |(/3-l))J 

i(/ 3 -l + a)J 


III. 


|*0.l) - 


1*1,3) 


IV. 


l*l,l)- 


1*0,2) 


V. 


l*l,l)- 


1*1,2) 


VI. 


l*l,l)- 


1*1,3) 


/3J 


VII. 


l*l,l)- 


1*2) 


i(/ 3 + l + a)J 


VIII. 


1*1,2) - 


1*0,2) 


(/ 2 - §a) J 
±(/ 3 + l-a)J 


IX. 


|*1.2) - 


1*1,3) 


X. 


1*1,2) — 


1*2) 


J 


XL 


|*0,2) — 


1*1,3) 


(-/2+!(/3 + l))J 

i(-/ 3 + l + a)J 


XII. 


1*1,3) - 


1*2) 



a This table uses the abbreviations fo = \J (2 — a s ) 2 + 3a 2 ,, 
/1 = Vl-a + a 2 , /2 = ^l-a/S + a 2 ^, /s = \/TTa T . 
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TABLE IV: Specific Heats" 



Spin System 


C/k B 


Dimer 


3(/3J) 2 e /3J /(3 + e' 3 ' 1 ) 2 


Symmetric Trimer 


|(/3J) 2 ei' 3J /(l + el' 3J ) 2 


Isosceles Trimer 


|(/3J) 2 (2(1 - a ) 2 e< 1+aa >" J + (2 + afe^+i^ 3 + 9a 2 ei Q ' 3J )/(2 + el"' 3 ' 1 + e (i+l«WJ) 2 


General Trimer 


±(/3J) 2 e^ 1+a ^ 3 (fZe^ 1+a ^ 3 + (4(1 + a s ) 2 + / 2 ) cosh(/ /3J/4) +4/o(l + a s )sinh(/ /3J/4))/ 
(l + e i( 1+a * )/3J cosh(/ /3J/4)) 2 


Tetrahedron 


18(/3J) 2 (lOe 2 "' 1 + 5e 3 ? 3 + ^ 3 )/($ + 9e 203 + 2e 3 '") 2 




" This table uses the abbreviation fo = \/ (2 — a s ) 2 + 3a 2 . 




TABLE V: Susceptibilities 2 


Spin System 


x/{<j^b) 2 


Dimer 


2/3/(3 + e l3J ) 


Symmetric Trimer 


i / 3(5 + ei^)/(l + ei' 3J ) 


Isosceles Trimer 


i/j(lO + ei a/3J + e ( 1 +3 Q '' 3J )/(2 + e i a " J + eP+s"" 57 ) 


General Trimer 


i/?(5 + e^ 1+a ^ 3 cosh(/ /3J/4)) /(l + e^+^/S-J cosh(/o/?J/4)) 


Tetrahedron 


2/?(5 + 3e 2 ' 3J )/(5 + 9 e 2 ' 3J + 2 e 3 ' 3J ) 


Rectangular Tetramer 


213(5 + e? 3 + e af33 + e^ +a ^)/(5 + 3e 03 + 3e a < 33 + 3e< 1+ ^ J + 2e< 1+ ^ J cosh(/ 1 /3J/2)) 


Linear Tetramer 


2/3(^5 + e 133 + 2e^ 1+a)f33 cosh(/ 3 /3J/2)) / 

(5 + 3e' 33 + 2e {1+ ? a)f33 cosh(/ 2 /3J) + 6e3 (1+a)/3J cosh(/ 3 /3J/2)) 



a This table uses the abbreviations fo = ^/(2 - a s ) 2 + 3a' 2 ,, fx = Vl - a + a 2 , f 2 = ^/l - a/2 + a 2 /4, / 3 = Vl + a 2 . 



